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APPLICATION OF RADIAL-EQUILIBRIUM CONDITION TO AXIAL-FLOW COMPRESSOR AND 

TURBINE DESIGN 

By Cbtjng-Hua Wtj and Lincoln Wolfenstein 


SUMMARY 

Basic general equations governing the three-dimensional com- 
pressible flow of gas through a compressor or turbine are given 
in terms of total enthalpy, entropy, and velocity components of 
the gas. Two methods of solution are obtained for the simplified, 

steady axially symmetric flow; one involves the use of a number 
of successive planes normal to the axis of the machine and short 
distances apart, and the other involves only three stations for 
a stage in which an appropriate radial-flow path is used. 
Methods of calculation for the limiting cases of zero and infinite 
blade aspect ratios and an approximate method of calculation 
for finite blade aspect ratio are also given. In these methods, 
the blade loading and the shape of the annular passage wall may 
be arbitrarily specified. 

The analysis shows that the radial motion of gas consists of 
a gradual, generally monotone component due to the taper in the 
passage wall, and an oscillatory component due to the radial 
variation of the specific mass flow at different stations along 
the axis of the machine specified in the design. The streamline 
is curved by this radial flow and a corresponding radial pressure 
gradient is required to maintain this curvature. The magnitude 
of this gradient is increased with high Mach number of gas 
flow and high aspect ratio of blade row. The conventional 
method of calculation, in which the effect of radial motion on 
the radial distribution of gas state is neglected, is found to be 
applicable only for the limiting case of zero aspect ratio. 

An analysis of the equations governing the flow shows that a 
designer is free to prescribe a reasonable radial variation of one 
of the velocity components or other thermodynamic properties 
of the gas at any station within the blade region. The various 
ways of using this degree of freedom and the different types of 
design obtained are discussed. Numerical computations are 
then made for two types of compressor and one type of turbine. 
The results indicate that, even in the case of nontapered passage 
walls, appreciable radial motion occurs and the corresponding 
effects are of significant magnitude and should be considered 
in design. 

INTRODUCTION 

The design of a compressor or a turbine (either of which 
is referred to hereinafter as “a turbomachine”) may be 
divided into two phases. The first phase concerns the type 
of design to be used, or the determination of the most desir- 
able possible variations of velocity and thermodynamic prop- 
erties of the gas in planes normal to the axis of the machine 


between successive blade rows. The second phase concerns 
the design of blades that will give the desired variations of 
velocity and other properties of gas in these planes. In the 
first phase, the condition of radial equilibirum (that is, the 
radial component of the equation of motion) must be used. 
The flow of gas in a turbomachine is curvilinear; it is curved 
not only by the whirling motion of gas, but also by the 
radial motion of the gas (reference 1). The equation of 
motion then specifies the radial pressure gradient required 
to provide the centripetal force to maintain the curved flow. 



(a) 



(a) Stream surface over four stages of multistage turbomachine. 

(b) Intersection of stream.surfaco with plane normal to axis. 

(c) Intersection of stream surface with axial plane. 

Figure 1. — Stream surface over four similar stages of multistage turbomachine and inter- 
section of stream surface with planes normal to and containing axis of machine. 

In figure 1(a), a curved stream surface over four similar 
stages of a multistage turbomachine is shown and figures 1 (b) 
and 1(c) show the intersections of this stream surface with 
planes normal to and containing the axis of the machine, 
respectively. The radial pressure gradient due to the 
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whirling motion of gas is always positive; whereas that due 
to the radial motion of gas may be either positive or negative, 
depending on whether the curvature caused by the motion 
is inward or outward from the axis of the machine at the 
point of consideration. Even when the radial motion 
involved is small, if the gas velocity is high and the blade 
aspect ratio is large, the radial pressure gradient due to 
the radial motion is of significant magnitude compared with 
that due to the whirling motion of gas and should be included 
in the design calculation. 

In the calculation of the state of gas in the normal planes 
far upstream and downstream of a single row of blades, 
where the radial motion is small, the pressure gradient is 
essentially due to the whirling motion alone. Experimental 
measurement checked well with the calculation when only 
the whirling motion was considered (references 1 to 3). 
For the general case of the gas in the normal planes between 
closely spaced successive blade rows, however, no satisfactory 
theory exists to calculate the magnitude of the radial dis- 
placement of the streamlines and its effect on the radial 
distribution of the state of the gas. A preliminary theoretical 
investigation of this problem conducted at the NACA 
Lewis laboratory was completed in April 1948 and is 
presented herein. 

In the analysis, the general equations governing the three- 
dimensional flow of gas in turbomachines are expressed 
in terms of total enthalpy, entropy, and velocity components 
of the gas. They are developed primarily for the case 
of steady axially symmetric flow corresponding to the limiting 
case of an infinite number of blades. Two numerical 
methods of solution are presented; one uses a number of 
successive stations through the turbomachine, the other 
uses only three stations for a stage in which an appropriate 
radial-flow path is employed. 

Methods of*^ solution for the limiting cases of very small 
and Very large blade aspect ratio- are then discussed. An 
approximate'^o^W of the radial displacement across a 
blade row having a finite aspectj-atio is given for the general 
case in which the whirling velocity 'of gas*is prescribed in 
design. 

The basic equations obtained are also used to investigate 
the maximum compatible number of radial variations of 
the velocity components and other thermodynamic prop- 
erties of gas that a designer is free to specify. It is found 
that the designer can specify only one such variation at each 
station along the axis of the machine within the blade 
region. Various ways of specifying this variation and 
the different types of design obtained are discussed. 

The methods developed are applied to two types of 
compressor and one type of turbine, in order to investigate 
the magnitude -of the radial motion and its effect on design 
calculations. 


FORMULATION OF EQUATIONS 

GENERAL BASIC EQUATIONS 

The three-dimensional compressible flow of gas through 
a turbomachine is governed by the following set of general 
basic equations (references 4 to 6) : 

From the principle of conservation of matter, the equation 
of continuity is 

i~+v.(pp)=o a) 


(Symbols used in this report are defined in appendix A.) 
The principle of conservation of momentum is expressed 
by the Navier-Stokes equation as 

P D j^ t = P F-Vp+pV 2 V + § V(V- P)+2[(Vp) ■ V]P+ 


(V M ) X(VXF) -f (V • V) (V M ) (2) 


where F is the external force exerted on unit mass of gas. 
The principle of conservation of energy may be written as 


Du 

Dt 

where u is related to T 


+P 

by 


Dip- 1 ) 

Dt 



Du__ DT 
Dt Cv Dt 


(3) 


(4) 


when conduction only is considered, Q is given by 

Q=p-^-(kVT) (5) 

and $ is the dissipation function given by 

$= M |2Vd(V-V)in+(VXT0 2 -2(P-V)(V.V)-| (V-P) 2 | (6) 

For the range of gas temperature and pressure usually en- 
countered in turbomachines, p, p, and T are accurately 
related by the following equation of state : 

p=R P T (7) 

Theoretically, the preceding seven equations, together 
with the given body force, known variations of c„, p, and k 
with temperature, and suitable boundary and initial condi- 
tions, completely determine the flow of gas through the 
turbomachine. It is found convenient in the present in- 
vestigation, however, to base the calculation on total 
enthalpy and entropy, which are defined by 


H=h+±V 2 ( 8 ) 

where 

h=u J rpp~ l (9) 

and 

Tds=du+pd(p - 1 ) (10) 
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By use of equations (8) to (10), the following forms of con- 
tinuity, motion, and energy equations are obtained (ap- 
pendix B): 


^ v +~xm^‘ T -m(S) = 0 (la) 


VH=F+TVs + VX (VXF)-^+ 


- V 2 V'+-V(V-F) 

PL 3 J 


+ 


^|2[(V M )-V]y+ (Vja)X(VXF) — | (V-F)(V M )| (2a) 

w^ Q+ t + Wt +r ( J+ p [ V27+ I v ( v ^')]+ 

^2l(V / x).V]T+(V M )X(Vxy)-| (V-F)(V M )J) (3a) 


p! 

Dt~T^ u p 


(3b) 


Equation (la) gives the continuity relation in terms of 
velocity, temperature, and entropy of gas. Equation (2a) 
relates the gradient of total enthalpy with body force, 
viscous forces, velocity, and other properties of the gas. 
This vector equation gives three scalar equations in three 
dimensions. Equation (3a) gives the rate of change of total 
enthalpy of gas along a streamline in terms of rate of heat 
additions, rate of work done by body and viscous forces, 
and so forth. Equation (3b) gives the rate of change 
of entropy along a streamline in terms of rate of heat con- 
duction and of dissipation of energy due to viscosity. 


STEADY AXIALLY SYMMETRIC FLOW 

The solution of the preceding general equations with a 
given set of suitable boundary and initial conditions is 
extremely difficult. Useful results may be obtained by 
considering, as first done by Lorenz in steam-turbine theory 
(references 7 and 8), the limiting case of an infinite number of 
infinitesimally thin blades. In this simplification, the force 
exerted on the gas by a blade element at any radius is con- 
sidered to be uniformly distributed over the stream sheet be- 
tween two neighboring blades at that radius, and is considered 
the body force F in the previous equations. Eor incompressible 
and frictionless flow, the value thus obtained gives an average 
value in the circumferential direction, provided the departure 
from the average value is small (reference 1). Because the 
number of blades is usually large, this simplification is con- 
sidered to be reasonable and is also used in the present in- 
vestigation. For steady inlet and exit conditions, all 
partial derivatives with respect to angular coordinate 9 
and time t are then equal to zero and the ^ate of gas is a func- 
tion of r and z only. 

The ideal case of a nonviseous gas will be considered first. 
In this case, there exist two more relations defining the 
problem. One is the fact that blade force is normal to the 


surface of the blade and, consequently, to the relative 
velocity of gas or the relative stream surface; that is, 

f-(v-U)== o (ii) 

or, referring to absolute cylindrical coordinates r, 6, z and 
the relative angular coordinate x, 

F r dr-\- rF g dx-\-F z dz=0 (11a) 

The other is the condition of integrability, of the. blade 
surface, 

F-(VxF)=<r (12) 


which in the case of axial symmetry reduces to (references 8 
and 9) 


J) 

dr 



(12a) 


From the general equations (la), (2a), (3a), and (3b), and 
equations (11) and (12a), the following equations are 
obtained for steady axially symmetric flow of nonviseous 
gas: (See appendix B.) 


1 d(rU r ) dUz 1 
r dr dz ‘7— 1 


(L^log. T+V^log, T)~ 




(13) 


, T , Ve d(rUa) 

dr T ' dr ‘ r dr ^ 2 \ dr dz J 


(14) 


0 = F»- 


1 ~ 

r J_ 


v d(rV 0 ) [ ^ b(rU 0 )1 


dr 


dz J 


(15) 



j) 

dr 



(12a) 


In the preceding equations, equation (13) is the continuity 
equation; equations (14), (15), and (16) are the three equa- 
tions of motion in the radial, circumferential, and axial 
directions, respectively. Equation (17) is considered to 
represent the energy equation and equation (18) to represent 
equation (11). In these equations, Q is now the heat trans- 
fer from the blade to the gas, uniformly distributed in the 
circumferential direction, as is the blade force F. These 
seven equations are considered seven independent equations 
that relate the eight unknown variables, which consist of 
three blade-force components, three velocity components, 
and H and s of the gas. The first three quantities determine 
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the shape of the blade and the last five quantities completely 
determine the state of the gas (all other thermodynamic 
properties of gas, such as p, p, and T, can be computed 
from them by using equations (7) to (10)). 

For compressors and turbines without blade cooling, the 
heat transfer between blade and gas is negligible; the en- 
tropy of gas is then constant along any streamline according 
to the energy equation (17). If the inlet air has a uniform 
value of entropy, the radial and axial derivatives of entropy 
in the preceding equations equal zero. 

In the case of real gas, the axially symmetric simplified 
forms of the viscous terms in equations (la), (2a), (3a), and 
(3b) can be obtained in a similar manner. These terms in 
the equations of motion may be neglected when compared 
with other terms in the same equation if the boundary layers 
along the passage walls are relatively thin. Because of the 
viscous shearing stresses in the gas adjacent to the blade, 
the force exerted by the blade on the gas is now slightly 
inclined from the direction normal to the relative velocity 
of gas and, consequently, equations ( 11 ) and ( 12 ) are not 
strictly true (the force components in the equations should 
be replaced by the direction cosines of the normal to the 
blade surface). Without using equation ( 11 ), however, 
equation (18) can be obtained from the equation of motion 
and the energy equation for steady flow with the assumption 
that the heat generated from the frictional work remains in 
each stream sheet (appendix B) and can therefore be con- 
sidered as representing the energy relation in the set of 
equations. The entropy increase along the streamline is then 
computed from a consideration of the actual compression or 
expansion process 

p=K P n (19) 


The formula (appendix B) 



In equation 

n may be directly obtained 
perature data. In a new design, 
the assumed poly tropic efficiency used in 
for uncooled blades: 

For compression, 


y ( 20 ) 

a given machine, 
ure and tem- 
btained from 
r gn calculations 



or 


1 7 — 1 
V Y 


1 - 


For expansion, 


( 21 a) 


may be adequate to account for the viscous effect in calcu- 
lating the pressure and density change along the streamline 
for the present problem. This correction is more important 
in the case of multistage compressors. For a viscous fluid, 
equation (20) therefore replaces equation (17) in the previous 
set of equations; and equations (14), (15), (16), (18), and 
(12a) are considered approximately true. With equation ( 12) , 
there are still seven independent equations defining the 
flow and the shape of the blade. 

METHODS OF SOLUTION 

The preceding section presents seven independent equa- 
tions relating the eight dependent variables V r , Ve, V z , H, s, 
F r , Fe, and F z , which define the flow of gas and the blade 
shape in the blade region. In the direct problem with a given 
machine, the shape of the blade section provides one more 
relation between Fe and F z , giving eight relations to deter- 
mine the variation of the eight quantities throughout the 
blade region. In the inverse problem, an appropriate desira- 
ble variation of any one quantity is prescribed within the 
blade region; the preceding seven equations then determine 
the variation of the remaining seven quantities throughout 
the blade region. No general solution of these equations 
seems possible, however, in either problem. Two numerical 
methods of solution are therefore suggested. In the first 
method, the preceding equations are applied to successive 
planes normal to the axis of the machine and short distances 
apart; this method is applicable to both direct and inverse 
problems. In the second method, a particular case is con- 
sidered, in which a simple appropriate radial-flow path is 
prescribed in the design. This method may also be used as 
a simple approximate solution in a direct problem in which 
the radial-flow path is approximately known. 


METHOD OP FINITE DIFFERENCE FOR SUCCESSIVE AXIAL STATIONS 


When two successive stations j and k a short distance apart 
within the blade are considered (fig. 2 ) and rV e is denoted 
by f, equations (14) to (16) may be written for each station 
as 




dr 


dr r 2 dr 


\ dr dz ) 


(14a) 


d_H 

dz 


o=fi-)(v.|+v.|) 


F,+ T~_+L¥.~Vr 


dz r 2 dz 


\ dr 


ayA 

dz ) 


(15a) 

(16a) 


The change in total enthalpy and entropy between the two 
stations is obtained from equations (18) and ( 20 ) (appendix B): 


a— 1 
n 

or 

Y— 1 

Y 



( 21 b) 




]+ £ Q dt (18a) 

tt Vn* 


n — 7 


(n 1 ) ( 7 — It 


log e 


V 2 

TT v ) 
tLj 


Because the change in s is usually small compared with the 
changes in H and V, the preceding method of determining s 


s t (r k ) — Sj(rj) =R 


( 20 a) 
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where (r) indicates that the gas properties at a particular 
station are a function of the radial position of the gas 
particle in that station. (It should be noted that because of 
radial motion, the radial position of a gas particle at any 
station k is different from its radial position at the previous 
station j.) 

Instead of integrating equation (13), the continuity 
relation between the two stations is readily obtained by 
equating the mass flow at the two stations: 


G k r k dr k = Gfj dr } (22) 

By expressing G in terms of H, V, and s, equation (22) 
becomes (appendix B) 

1 Si 1 >£ 

R r } dr } 
(22a) 

Equations (12a), (14a), (15a), (16a), (18a), (20a), and 
(22a) are now seven independent equations relating the gas 
properties and the blade forces at the two stations in the 
blade region. In these equations, the heat transfer is 
negligible in an ordinary turbomachine and can be estimated 


V, 




e r k dr k ~V 2t j ( Hj — 




in the case of cooled turbine blades; the temperature T is a 
known function of H, V, and s; tv is given; and r k is obtained 
from rj and V rJ . Hence, there are only eight unknowns in H, 
s, V r , Ve, V„ Fr, Fg, and F z at the second station k. In the 
direct problem, the blade shape gives one more relation 
between Fg and F z \ whereas in the inverse problem, one 
suitable relation among the eight unknowns is specified by 
the designer. In either case, the unknowns at station k can 
be obtained from the known values at station,?', the passage- 
wall shape, and the preceding relations. 

In the free space between two blade rows, or in the space 
upstream of the first blade row and downstream of the last 
blade row, the force terms drop out of the equations, which 
results in equation (22 a) and the following equations 
(neglecting friction and heat transfer between gas and 
passage wall) : 


5s , /ay 2 bV r \ 

dr ()rr 2 dr ' 2 \ dr 5s/ 


0=V|+F, 


5f_£r 
5s Dt 


5I?_ T 5s , ! bl_ v /5F, bV T \ 

5s 1 5s ' r 2 ds r \ 5r 5 z ) 


(14b) 

(15b) 

(16b) 

(20b) 


Equation (22a) and these four independent equations, 
together with the given passage-wall shape, completely 
determine the variations of the five independent quantities 
V T , Vt, V t , H, and s outside the blade region. The solution 
of the problem over the entire region inside and outside the 
blade region, using this step-by-step method, varies with the 
type of design and the condition given or prescribed. In 
any case, the computation would be quite laborious. 

In order to obtain an over-all pictur^pf the-radial flow in a 
turbomachine and its effect on^design consideration in a 
simpler way, the following method considering the problem 
of a particular case is given: 


METHOD OF PRESCRIBED RADIAL-FLOW PATH 

In a turbomachine, the radial motion of the gas is caused 
by three factors: 

(1) Tapering of the annular passage either at the inner or 
outer wall gives the flow a radial displacement across the 
stage, which is, of course, greatest in the immediate neighbor- 
hood of the tapered surface. 

(2) Even with a nontapered passage, a radial displacement 
across the stage may be necessary because of a variation 
in the distribution of specific mass flow over the blade height 
across the stage. 

(3) Even if no radial displacement occurs across the stage 
(that is, the same particle occupies the same radial position 
at the first station of each successive stage), there will, in 
general, be radial displacement of flow within the stage. 
This radial flow will then be oscillatory in nature, a radial 
displacement in the rotor being followed by an equal and 
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opposite radial displacement in the stator. This radial flow 
arises because of the difference between the radial variation 
of the specific mass flow within the stage and that at the 
entrance and exit stations of the stage. (This radial displace- 
ment can only be avoided by specifying zero or the same radial 
variation of specific mass flow at all stations of the stage 
in the design.) 

In general, the radial flow of gas therefore consists of a 
gradual, generally monotone, radial motion due to factors (1) 
and (2), with an oscillatory motion of period equal to the 
stage length due to factor (3) superimposed on it. The 
radial flow caused by these three factors will be similar to 
that shown in figure 1. The effect of the radial motion 
on the calculations arises chiefly through the term bV T /dz in 
the radial-equilibrium equation (14a). This term is expected 
to be significant mainly because of the oscillatory motion, 
which may require significant changes in V T within a single 
row of blades. The case of oscillatory motion within a stage 
with no over-all radial displacement across the stage will 
therefore be considered first. That is, the gas-passage wall 
is nontapered and the radial distribution of gas properties 
at the entrance and exit stations of the stage is the same. 

Because there is no blade force acting on the gas, the gas 
flowing through the gap between two blades is under a 
nearly constant pressure gradient and consequently tends 
to move with the same curvature it acquires while leaving 
the first blade. For nontapered passages, the maximum or 
minimum point of the radial-flow path is likely to be some- 
where near the middle of the gap. (The intersecting curve 
of a stream surface with an axial plane is herein referred 
to as “the radial-flow path.” Because of axial symmetry, 
the radial-flow path is the same in any axial plane.) The 
stations between blade rows are most conveniently chosen 
at these points. The stations in front of the rotor, between 
the rotor and the stator, and behind the stator are denoted 
by subscripts 1, 2, and 3, respectively. (See fig. 3(a).) If 
r 0 and L represent the mean radial distance of the flow path 
and the axial length of the blade row, respectively, then the 
radial distance of the gas particle at position z is given by 


r~r 0 


r 2 —r i 


2 



(23) 


and at stations 1, 2, and 3, z/L=0, 1, 2 


/(0)=/(2) = i,/(i) =— n 

./' ( 0 ) =/'(!)=/' ( 2 ) =0 ) 


(24) 


w^-re/ is a function giving the form of the radial-flow path 
and the prime indicates differentiation with respect to z/L. 
It follows that 


1 2 3 



1 2 
U-i— *| 



(a) Nontapered passage. 

(b) Tapered passage. 

Figure 3. — Stations between blade rows. 


At station 1, 2=0, 


and 



V tA = 0 


r 2 —r x 
2 L 2 




(26) 


inasmuch as bVJdz is practically zero in passing through 
the gap. Similarly, at station 2, z—L, 


and 


V,, 2 = 0 


/dFA r 2 —r x 

\ dz J 2 2 L 2 


F, 2 /" (1) 


(27) 


Because j"(z/L) determines bVrjbz or the effect of radial 
motion on the radial-equilibrium condition, it is desirable 
that it vary continuously; this condition together with 
those of equations (24) suggests 

y(|)= c °sx(|) 




r 2 —rj 
2 L 


V, 




(25) 


Then 


/"( 0) = — * 2 , /"( 1)=7T 2 
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and equations (23) and (25) to (27) become, respectively, 


r ~ r o— 2 2 1 cos 7T ( z ) 

(23a) 

v,=*^v,sin r (£) 

(25a) 

/dFA 2 r 2 — r, T/ 

fdzjr' 2 u 

(26a) 

(Wr _,r w v ‘ * 

(27a) 


For the sinusoidal form of f{z/L), the maximum absolute 
value of f"(z/L) occurs at 2 = 0 , L, and 2 L, and is equal 
to t 2 . Even if f"(z/L) is assumed constant between 2=0 
and z—L/2, thus minimizing the maximum absolute value of 
f" in the interval, the absolute value off" equals 8. This 
assumption, however, necessitates a discontinuity in j" at 
z—L/2. The values of 7 r 2 for the absolute values of/"(0) 
and /"(l) can therefore be considered as small as is likely. 
The smooth variation of f'iz/L ) and the minimization of 
the absolute value of f"{z/L) at the stations make it desir- 
able to employ equation (23a) in the design. This simple 
radial-flow path will also give a good approximate answer 
to most designs in which the distribution of blade loading 
and radial blade forcfe is not too ununiform in the axial 
direction. In such cases, this simple sinusoidal curve is 
believed to represent the major harmonic of the actual radial- 
flow path and the principal effect of the radial motion may 
be obtained through the use of this simple curve. 

The radial-equilibrium equation (14b) may be written in 
terms of (r 2 — rO by use of equations (26a) and (27a) 


(IHi rrt 

W i ~~ Ll 


du 

dr, 'r, 2 dr, 


+v t ,. 


dV,, t 
dr , 


+ (-DV 2 


r»—r 1 
2L 2 


(14c) 


where 


i= 1, 2, 3 


When r, is replaced by the dimensionless variable 



equation (14c) becomes 


dH, rp ds , 1 1 r i d / \ 1 

drf~ * dr,r , 2 r t —r h dr / \r t — rj^ 


V, 


dr 1 


'+(-!)* w (rf-r^AW,., 2 (14d) 


This form of the radial-equilibrium equation is seen to con- 
tain a term directly proportional to the radial displacement, 
to the square of the axial velocity, and to the square of the 
blade-row aspect ratio. If the blade-row aspect ratio is 
large or the axial velocity is high, the effect of radial motion 
may be large even though only a small amount of radial dis- 
placement occurs across the blade row. 


This method is readily extended to the case where an over- 
all radial displacement occurs across the stage due to tapering 
of the passage or due to variation in the design from stage 
to stage. In figure 3(b), the radial position of a gas particle 
originally at r, in station 1 is at r 3 in station 3. For the 
oscillatory motion required within the stage, r 2 is not gener- 
ally equal to ^ ( 7 * 1 + 73 ). For the same reason stated in the 

previous case, it is desirable to have the radial-flow path 
consisting of a sinusoidal curve superimposed on the line 
passing through ( 2 1 , ri) and ( 23 , r 3 ); that is, 


— +^1 +K^X 1 - C03 7) (23b) 

Then 

and 

(26b) 


inasmuch as dVJdz is practically zero in passing through 
the gap. With this value of bV r /c)z, the radial-equilibrium 
equation (14b) becomes 


dHj rp ds i . i dfj . 

dr , ~ 1 dr,ri 2 dr, z,i 


dV 2 ,i 

dr , 


+ (-l) 




n+rs 

2 


) 


V, 


(14e) 


This equation is similar to equation (14c). (If r^=r u it 
reduces to equation (14c).) A similar equation in dimen- 
sionless r / can also be obtained for this case by dividing r, 
by (r t — r h )i. 

With this method of prescribed radial-flow path, if it is 
only required to find the velocity distribution in stations 
between successive blade rows (to provide data for the design 
or setting of blades), the distribution can be obtained by 
considering only these planes without making any computa- 
tion in the blade regions. For example, suppose that all quan- 
tities are known at the inlet station of the stage (station 1 ), 
then for station 2 behind the rotor blade row the follow- 
ing relations exist: When the untapered passage walls are 
considered, equations (14c), (18a), (20b), and (22a) give 


n — 7 


and 


V,. a (ffa -5 vv) 


1 $2 
-i ~i ~R , 
e To di\ 


2 T - \T 2 

~2JF Vz ’ 2 

(14f) 

t- J) Q dt 

(18b) 

fla-i W 

( 20 c) 

Hy—Vf 

\ 

$1_ 

R r 1 dr. 


( 22 b) 


873599—50 2 
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If the radial variation of one quantity is known (such as 
&r 2 /&r 2 = 0 , dVg i2 /dr 2 — 0, or dj8 2 /dr 2 =0) among the five un- 
knowns at station 2 (consisting of f 2 , V Zi2 , II 2 , s 2 , and ^(ri)), 
all the other quantities can be determined. The procedure 
will vary with different types of design. In most cases, the 
computation may be started with an appropriate value of r 2 
as a function of r x obtained by an approximate solution. 
With the variation of one of the four quantities H 2 , s 2 , f 2 , and 
F 2i2 given, the remaining three quantities are computed from 
equations (14f), (18b), and (20c). These values are then 
inserted in equation (22b) to determine if the continuity 
relation is satisfied. If it is not, the values must be adjusted 
until the continuity relation is satisfied. An alternative pro- 
cedure is to insert these values in equation (C 5), which is 
derived from equation (22b), to obtain new values of r 2 as 
a function of r x , and the entire process is repeated until the 
desired accuracy is reached. In a tapered passage, both 
stations 2 and 3 must be computed at the same time. An 
application of this method to a nonvortex-type compressor 
stage will be given in the section NUMERICAL EXAMPLES 
FOR TYPICAL DESIGNS. 

LIMITING CASES AND APPROXIMATE SOLUTIONS 

Simplified -radial -equilibrium approximation. — In this 
commonly used approximation, the gas is assumed to flow 
on cylindrical surfaces for nontapered passage walls; that is, 


V T =0 

r x =r 2 =r z 


(28) 


At these stations, the radial-equilibrium equation (14b) 
reduces to 


(ill j m dSj . £ j i ■ y y dV Z< l 

dfi i dri r? dri z,i dvi 


(14g) • 


where i= 1, 2, or 3. 


For tapered passage walls, a certain simple relation among 
r u r 2 , and r 3 is assumed, but the term containing dV r j dz in 
the radial-equilibrium equation is still neglected (that is, 
equation (14g) is used instead of equation (14b)). 

With this simplifying assumption, the gas state at station 
2 or 3 can be computed much more easily from the gas state 
at station 1 and the one condition specified at stations 2 
and 3. The continuity equation (22b) for individual stream 
sheets, however, is now discarded because of the assumed 
relation among r u r 2 , and r 3 and is replaced by the following 
continuity relation for the entire annular area: 



(22c) 


Thus, equations (14g), (18b), and (20c) are used in this 
calculation, with equation (22c) used as a check on total 
mass flow. 

When dsijdri is negligible, the following equations for the 
radial variations of pressure and density can be obtained 
from equations (14g), (B4), and (B8): 


Pt,i y „ ~y dPi^Y 9 J dVr ’j 
Pt.i Pt dr t r t r,i dri 

P t,i 7_ 2 dpi Ff),i tt- d V r ,i 
7 Pt,i y Pt dr~ r t rA dr t 


(29) 

(30) 


In these two equations, the last term is very small compared 
with the next-to-last term and may therefore be neglected. 
Under the present assumption, this term becomes zero for 
a nontapered passage and the resulting equation may be 
more directly obtained, as is usually done, by taking the 
approximation involved in the use of 

1 dpj Ve, i 
Pi dri Ti 


for the equation of motion in the radial direction in place of 


1 dpi 
pi dri ' 


V, 


Ti 


y f<>Z.A_Tr d JYt 
■ *• 4 \dz )i '■ 1 dTi 


Equations derived from this basis of calculation in dimen- 
sionless forms are given in appendix D for several types of 
design. 

Limiting case of zero aspect ratio. — Two limiting cases will 

now be discussed for which the evaluation of the term ^ 

can be avoided. If the blade row has an axial length suffi- 
ciently great relative to the radial length (that is, if the 

blade-row aspect ratio is sufficiently small), the term ^ Dz\ 

will be negligible in spite of any radial displacement across 
the blade row. This extreme situation is designated the 
zero-aspect-ratio case and differs from the simplified -radial- 
equilibrium approximation in that the radial displacement 
across the blade row is properly determined and its effect 
on the state of gas is included in the calculation. The 
continuity equation (22b) for individual stream sheets is 
therefore to be satisfied in addition to equation (22c), but 
equation (14g) is still used in place of (14c). In the case 
of gas-passage walls having no taper or slight taper, the 
difference between the two cases is small; a successive- 
approximation procedure starting with the result of the 
simplified-radial-equilibrium calculation can therefore be 
used. This procedure may be outlined as follows: 

1. With the given values at station 1 and one prescribed 
condition at station 2, use the simplified-radial-equilibrium 
equations to find f 2 , V z 2 , H 2 , and s at station 2 as functions 
of r x ; then compute G 2 (r l ). 

2. By using the value of G 2 (r x ) obtained from step i, find 
r 2 (r 1 ) from equation (C5). 

3. Substitute this value of r 2 (r t ) into equations (14g), 
(18b), and (20c) to obtain a second solution for f 2 , V Zi2 , FL, 
and # 2 as functions of r u 

4. Repeat steps 2 and 3 if necessary, using the value of 
G 2 (ti) obtained from step 3. 

In the case where there is considerable taper at the passage 
walls, it is better to assume r 2 {r x ) according to the taper to 
start the calculation rather than to use steps 1 and 2. 
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Limiting case of infinite aspect ratio. — The other limiting 
case corresponds to a blade row with axial length negligible 
as compared with radial length, and is designated the infinite- 
aspect-ratio case. The negligible axial length does not 
provide space for any appreciable radial displacement, 
hence r 2 may be taken as equal to r u or 

G 2 =G X (31) 

and for a tapered passage, either the preceding equation or a 
relation similar to the one that follows may be used: 


functions of r u an approximate solution for the radial dis- 
placement across the blade row can be determined in the 
following manner: 

First take as two separate functions 

A e (ri) the function (r 2 —r x ) of r x satisfying the radial- 
equilibrium and total-enthalpy-change equations for 
a given distribution of the other variables 

A c (r x ) the function (r 2 —r x ) of r t satisfying continuity equa- 
tion (22) for a given distribution of other variables 


G 2 —G x 


r x ,t 2 —r x . h 2 

r 2 . t 2 —r x . h 2 


(32) 


Either of these two equations now takes the place of the 
continuity equation (22b). 

Because the change in axial length for a very small change 
in V r is also very small, (dTft/dzft does not vanish. Although 
its absolute value does not affect the radial motion because 
of the negligible blade-row axial length, the relative value of 
dV r /dz ahead of and behind a blade row is needed to 
determine completely the distribution of gas properties at 
these stations. If the loading of the blade in the axial direc- 
tion is relatively uniform or the blade is designed to give a 
sinusoidal radial-flow path, the curvatures of the radial- 
flow path at the two stations can be considered equal in 
magnitude and opposite in sense. Then, 



(33) 


In order to combine this relation with equation (14b) in a 
simple manner, it may be assumed that 


It is assumed in this method that the radial gradients in 
V z and f depend primarily on the magnitude of the radial 
displacement ( r 2 —r x ) and not on its exact distribution. 

Accordingly, 

A eGO^Ve g(r x ) (35) 

where y e is the maximum value of A e and g{r x ) is a plausible 
form for the distribution of A e satisfying the boundary 
conditions : 

g(ri,h)=g(r Xwt )=0 ) 

for [ (36) 

g' (?i.n) = 0 , g{r x ,^) = 1 ) 

If A<. is calculated for several values of A e , it is possible to 
plot y c , the maximum value of A c , against y e ■ A fairly good 
approximate solution might be expected to correspond to 
the point y c =y e . This process can be further refined by 
varying g(r x ) from the function originally assumed in the 
direction of the calculated function A Jy e . 

By the use of this procedure, the following approximate 
value for the magnitude of radial displacement across a 
blade row is obtained (appendix F) : 


V itX 



(34) 


r 2 r i === 


(r 2 —r x ) 
1 +A 2 


s 


(37) 


Combining equation (34) with equation (14b) at stations 1 
and 2 yields 

(IH X . (Ill 2 rp ds 1 . rp ds 2 . x rfft . ft d% 2 I 

dr x ' dr 2 ~ 1 dr x 2 dr 2 r x dr x r 2 dr 2 


This value can be used as a starting value for exact calcula- 
tions or may be used as the final value for approximate 
calculations. 

APPLICATION TO DESIGN 

DEGREE OF FREEDOM IN DESIGN 


V, 


dV ZiX | T r dV z,2 

K^ +V - 2 dr 2 


(34a) 


For a typical stage of a given design, equations (14b), 
(18b), and (20c), with either equation (31) or (32) and equa- 
tion (33) or (34) will completely determine the variation of 
gas properties at the two stations. 

In appendix E, formulas are given in dimensionless forms 
for two common types of design in order to calculate the 
variations of gas properties for the two preceding limiting 
cases. The results so obtained will give the limits of the 
variation of the gas properties along the blade height. If 
the difference is large, it is worthwhile to make the calcula- 
tion for the given blade-row aspect ratio. 

Approximate solution for finite aspect ratio. — For the 
general cases where ft and ft are prescribed in design as 


In the preceding analysis it was shown that within the 
blade region there are only seven independent equations 
relating the eight dependent variables that determine the 
state of gas and the blade shape; whereas in the space outside 
the blade, there are five independent equations that deter- 
mine the five dependent variables determining the state of 
gas. In the inverse or design problem, the designer therefore 
has one and only one degree of freedom for prescribing a 
reasonable radial variation of one single quantity in all blade 
regions, of a different quantity in different blade regions, or 
a single relation between several quantities. In addition, he 
is free to specify the taper in the passage wall, the position 
where the radial element of the blade is set, and a suitable 
condition of gas at stations far ahead of and behind the 
machine as boundary conditions, such as a uniform state of 
gas entering the machine. 
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In current design practice, design computation is made 
only in the planes between successive blade rows in which 
a certain desirable variation of one gas property can be 
specified. The blades are then either selected from cascade 
data or theoretically designed to achieve the calculated gas 
states in these stations. For this reason, the following dis- 
cussion of different ways of specifying these variations will 
be centered in these planes: For simplicity, only the adiabatic 
case is considered. Under this condition, if equation (18a) 
is applied to the three successive stations of a stage and is 
differentiated with respect to r x , 

dHx , . . 

dr i ' “ \^2 dr i dri ) dr 2 dr x dr z di\ ^ 


This condition is commonly used in turbines and compres- 
sors. Ignoring radial motion, in addition to this condition, 
constant total enthalpy and constant axial velocity over the 
blade height can be obtained. Considering radial motion, 
only one of these two additional conditions can be obtained 
in conjunction with equation (41). (See section NUMER- 
ICAL EXAMPLES FOR TYPICAL DESIGNS.) 

(4) Symmetrical velocity diagram. If V z j=V Zy2 and 
r l = r 2 =r, the symmetrical velocity diagram gives 

U<M + Ve t %=(ar (42) 

or 

£i + £ 2 = &r 2 (42a) 

Differentiating with respect to r yields 


A few ways of taking up the degrees of freedom at these 
stations between successive blade rows are listed in the 
following paragraphs: 

(1) Constant work per unit mass of gas flow over the blade 
height. This condition is usually specified in the design of a 
turbomachine. It relates £ behind the rotor to its value 
ahead of the rotor by 


£2(7*2) — £1(7*1) 

d £ 2 dr 2 d£ t 

dr 2 dr x ~~ dr x 


(39) 

(39a) 


where r Xft 8 t U ut is equal to (£ 2 — £ 1 ) at the blade tip and is 
also equal to (£ 2 — £ 1 ) at other radii. 

Constant work over the blade height gives constant total- 
enthalpy change over the blade height. If the velocity at 
the exit of a stage is equal to that at the entrance, this 
condition also gives constant static-enthalpy change over 
the blade height. 

Under the condition of constant work, equation (38) 
reduces to 

dH x _dJl 2 dr 2 _dH 3 dr% 

dr x dr 2 dr x dr 3 dr x ■ a ' 


(2) Constant total enthalpy over the blade height: 



(40) 


f+#= 2 “- ■ < 42b > 

If r x ^r 2 or V ZiX ^V Zf2 , the symmetrical velocity diagrams 
may be defined by 


Then, 


£i(ri) +£2(7*2) =wri 2 


d£± , dfr 

dr x dr x 


= 2w r x 


(431 

(43a) 


With the use of the symmetrical velocity diagram, the aero- 
dynamic limitations of gas flow through the rotor and the 
stator are reached at about the same time. Reference 10 
shows that the blade-profile loss is a minimum with the 
symmetrical velocity diagram if the lift-drag ratio is con- 
stant. For incompressible flow, the change in static pres- 
sure or enthalpy is also the same in passing through the rotor 
or the stator, and the stage is therefore often referred to as 
the “50-percent reaction stage.” 

(5) Wheel-type distribution of tangential velocity: 


V e , t =K i r i 

(44) 

or 


— — 2 K r 
dr t ~ Ki i 

(44a) 

(6) Constant tangential velocity: 


V 0A =K t 

(45) 


This condition usually applies to the first stage of a com- 
pressor and will hold for all succeeding stages if constant 
work per unit mass over the blade height is employed. If a 
nonzero value of dH x jdr x is desired, an initial preparatory 
stage must be specially designed to obtain this value. In 
the last stage, however, it is usually desirable that dH/dr be 
nearly zero. 

(3) Free- vortex-type distribution of tangential velocity: 



(41) 


dXi 

dr t 


=K t 


(7) Constant axial velocity over blade height: 


dr 


-0 


(45a) 


(46) 


At a very low speed of gas flow with no change in density, 
the specific mass flow will also be constant over the blade 
height ; there will therefore be no radial flow across the blade 
row and equation (14b) reduces to (with the entropy varia- 
tion neglected) 


dtf_£ d£ 

dr r 2 dr 


(47) 


or 


fi=K t 


(41a) 
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And when 


dH 

dr 


=0 



The equivalence of equations (46) and (47) breaks down, 
however, for current aircraft applications, where the speed 
of gas flow is high. 

If equation (47) is substituted into the radial-equilibrium 
equation (14b) with entropy variation neglected, the follow- 
ing relation is obtained: 


(9) Relative Mach number. For high performance, a cer- 
tain variation of relative Mach number consistent with the 
radial variation of solidity and thickness of blade may be 
specified in design. Then for the rotor, 


V^+iVe, i-urO 5 


(v-1) [h,-| (JV+f,., s )] 

For the stator, 

(7-1) (F 2 /+TV)] 


--M i 2 =<p 1 (r l ) (50) 


-M% — <p 2 (t 2 ) (51) 


d V z _dVr 
dr dz 


(47a) 


The left side of equation (47a) is the tangential component 
of fluid rotation VXF; thus equation (47) is a condition 
for potential flow in the free space between blade rows. 

If it is desired to take into account in design the effect 
of the boundary layers at the inner and outer walls of the 
gas passage, instead of equation (46) an appropriate axial- 
velocity variation close to the actual one may be prescribed 
in design: 

(46a) 


(10) Untwisted rotor blades. For simplicity in fabrica- 
tion, especially for a cooled turbine, untwisted rotor blades 
may be used. Inasmuch as the flow angle is only slightly 
different from the blade angle, the following relation may be 
used in design: 


fa 

dry 


=0 



(52) 


(11) Blades with all elements radial. For high-speed rotor 
blades, in order to reduce centrifugal stress it may be desir- 
able to have all blade elements radial. Then, 


(8) Constant specific mass flow over blade height. In 
order to avoid radial movement across the blade row in com- 
pressible flow, it has been suggested (for example, reference 11) 
that constant axial velocity be replaced by constant specific 
mass flow: 


dG\ dGz dG 3 

dr dr dr ~ 


(48) 


F T = 0 


By using this relation, equation (12a) reduces to 


or 


~=Kr 
Fo ^ 

tan f3=Kr 


(53) 


(53a) 


Radial displacement can also be prevented by the use of two 
conditions instead of three: 


dG\ dGi dG% 

dr dr dr 


(49) 


For designs using either of these two conditions, the simplified- 
radial-equilibrium calculation is more correct. Designs 
employing no radial flow have the advantage that the cal- 
culation does not involve any radial displacement across 
the blade row and that the two-dimensional-cascade data 
can be directly applied. The final equations derived from 
these conditions (equations (48) or (49)), however, are diffi- 
cult to solve and the conditions are incompatible with 
tapered passage m a multistage turbomachine. 


where K is a function of z. 

In multistage machines, similar variation in either tan- 
gential velocity, axial velocity, or specific mass flow may be 
specified at the. similar stations of each stage: 


#1 

dij 

(54) 

dr x 

dr 3 dry 

dV 2A 

dry 

_dV t , 3 dr 3 
dr 3 dr x 

(55) 

dGy 

dr x 

dG 3 dr s 
dr 3 dry 

(56) 


Stages of multistage machines designed for similar variations 
of gas properties from stage to stage arc termed “typical 
stages.” 
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Types of design. — A large number of different types of 
design may be obtained by different combinations of those 
conditions specified in equations (39) to (56). These 
designs may be conveniently divided into two groups. In 
the first group, the condition of constant work at all radii is 
specified in the design. That is, equation (39) is employed, 
which gives : 

d£ x _ d £ 2 

dr ~ dr i 

and 

dH x _JLH 2= dH z 

dr x , dr i dr x 



In cases where the symmetrical velocity diagram is also 
specified, by using equation (43a), equation (57) reduces to 


and 


dr x dr x 1 


dH x _dH 2 dH, 
dr x dr x dr x 


(58) 


In the second group, the condition of constant work is not 
specified. 

The following tables present a few types of design in each 
of the two groups for multistage turbomachines consisting of 
a number of similar stages. The manner in which the 
degree of freedom is used up at each station of a typical 
stage and the known characteristics of each type are given. 
The typical stage is considered to consist of a rotor followed 
by a stator. 

These conditions specified for the typical stage, as given in 
the table, completely determine the flow over all the stages. 
The flow in the inlet guide vane, if required, is to be deter- 
mined by the given condition at the inlet to the machine 
and the condition specified at station 1. Similarly, the flow 
in the last stator is to be determined by the condition speci- 
fied behind the last rotor and the given condition at the exit 
of the machine. 

It may be desirable for certain applications to use different 
types of design in a multistage unit. Those designs can be 
obtained by using only the relations of different types at 
stations 1 and 2 in the table. 



1. Free vortex 


rfr, 

Constant work 


dr, dr 2 dr 3 

dH l dH. 2 dH 3 
dr, dr, dr. 


V 9 —KjT; If —=0 at inlet to mac! 

Constant axial velocity over blade height for incompressible flow. dH 

Small radial gradient in axial velocity for compressible flow. — =0 at a11 sttttions - 


2. Symmetrical velocity 
diagram 


3. Wheel-type tangen- 
tial velocity in front 
of rotor 


Symmetrical ve- 
locity diagram 


Constant work 


Constant work 


i df* 2 d{ 3 


dHl dH i iH 3 
dr, dr, dr. 


<ff 1 #2 #3 , Kr 

dr, dr, dr, 1 

dHi dH 2 dH 3 
dr t dr l dr. 


„ "’’l S,U u ,r ut 

2 -“7T 

i <ar i ■ ** u u t r i. A r i 

^) v 2 


Combination of wheel-type and If — -=0 at inlet to machine, 
vortex-type tangential velocities. “ r 


Large negative radial gradient of axial velocity at all stations. 
F 9 ,, = Ar, 

Vt.t -(jfr,+a t cr, i( ^r ! ) ~ 

Large negative radial gradient of axial velocity at all stations. 


dH 

— =0 at all stations. 

dr 


-=0 at all stations. 


Constant work 


4. Same variation in 
axial velocity 


5. Same variation in 
specific mass flow 


dV zA 

dV, t2 




dr. 

dT \ 

dH x 

dH 2 

JJh 

dr7 


~ dr i 

dv i.l 

dV Zt 3 




dT \ 

dr x 

dV t 

.1 dV *,2 < 



dr. 

= dr 2 

Constant work 

dr, dr. 


da. 

, ja 2 


* 

= dr i 

dHi 

3 

JJh 

dr^ 

“dri 

dr x 

da. 

!_d0 3 




si 

"dr, 

da x 

d0 2 

dQ z 

dr. 

dr x dr t 


„ . ... _ ... . . . .. If —=0 at inlet to machine, 

For incompressible flow, this type requires no radial flow across dr 
blade rows and is equivalent to first type in group. For com- 
pressible flow, small radial gradient exists in f. °H 

H -—=0 at all stations. 


If — = 0 at inlet to machine, 

fo radial flow across rotor and stator blades for nontapered or 
passage. 

« — =0 at all stations, 
dr 
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Type 

Station 

Conditions specified 
at 3 stations of any 
stage 

Characteristics of flow 




dG, 












Constant specific mass flow 


dG, 

dG, 

dC 2 




dr s 

dr. 

dr 2 dr 3 

Ko .adial flow acioss all blades m uontapered passage. 



3 

d0 3 

dr 3 






1 

dV e , 
_-U _ 0 

dr, 




2. 

Constant axial velocity 

2 

dV. 2 
dr , 

dV z 

dr, 

, dV !i2 dV z-3 
dr 2 dr 3 

No radial flow across all blades in uontapered passage for incompressible flow. 



3 

o 

l-S 

"53 1 






1 

dr, 




3. 

Untwisted rotor blade 

2 

p-o 

dr 2 

d&i d & ^ d 0 3 

dr t dr 2 rfr 3 





dr 3 
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NUMERICAL EXAMPLES FOR TYPICAL DESIGNS 

The methods of calculation previously outlined are applied 
to the typical stages of compressors of types 1 and 2 of group I, 
as given in the table. The inlet total enthalpy is assumed 
uniform with respect to radius; and with work exchange 
with rotor uniform along the radius, total enthalpy is con- 
stant with respect to radius in all stations. The calculation 
is rendered dimensionless by expressing all velocities in terms 
of U t , total enthalpy in terms of U, 2 , and r in terms of r t . 
Because the main purpose of the calculation is to deter- 
mine the magnitude of the oscillatory radial motion and its 
effect on the radial distribution of gas properties, a non- 
tapered passage wall is used. Heat transfer is assumed to 
be zero in the calculation and the entropy is assumed to be 
constant at each station. The change of entropy across 
the blades at all radii is assumed equal to that obtained from 
the polytropic efficiency assumed at the mean radius. This 
calculation does not take into account the boundary layers 
at the rotor drum and the outer casing, and consequently 
is good only for the main portion of gas flowing between 
them. This restriction can be removed if more data on the 
variation of n with radius are available. 

In the comparison of different blade-row aspect ratios in 
each design, in addition to the same aerodynamic limitations, 
the same axial velocity at the mean radius is used. The 
comparison between different cases will be slightly different 
if another basis of comparison is used. 

Symmetrical-velocity-diagram and constant-total-enthalpy 
compressor. — Because the difference between zero- and 
infinite-aspect-ratio cases is found to be large in this design, 
two calculations are made for a blade-row aspect ratio of 2 ; 


one calculation is based on a prescribed sinusoidal radial- 
flow path, the other is based on the approximate solution of 
equation (37). The equations based on prescribed sinu- 
soidal radial-flow path are given in appendix G. The follow- 
ing design constants are used for all cases : 


Hub-tip ratio 0. 6 

Limiting Mach number relative to rotor blade 0. 8 

Limiting value of — 0. 7 

V 8,1 

Palytropic efficiency at mean radius 0. 9 

V'.uJUt - 0.772 


The limiting value of ^0 ■* _ ¥?, j j s b ase( f on a formula given 

V z.l 

by Howell (reference 12). The last value results from the 
use of in the simplified-radial-equilibrium 

calculation, and is used for all cases. The results of the 
calculation are shown in figure 4. 

The distribution of specific mass flow ahead of and behind the 
rotor for the different cases considered is shown in figure 4(a). 
It may be seen that in all cases except the infinite- 
aspect-ratio case, the specific mass flow G/G t increases to- 
ward the hub faster behind the rotor than ahead of the rotor; 
that is, passing through the rotor, the gas moves toward the 
axis of the machine. The magnitude of this displacement is 
obtained from the continuity equation (Co) and is shown in 
figure 4(b). In the simplified-radial-equilibriuin calculation, 
it is assumed that there is no radial motion, but when the 
distribution of specific mass flow is substituted in the con- 
tinuity equation (C5), quite large radial displacement across 
the blade is obtained. This kind of calculation is therefore 
inconsistent. In other calculations, the distributions of gas 
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properties are calculated from assumed radial displacements 
that are to be checked with the displacements required from 
the continuity relation with these distributions, and are 
therefore consistent in themselves. The radial displacement 
used in the approximate calculation for A=2 is obtained by 
the approximate formula, equation (37), and is about 25 per- 
cent lower than the value obtained from using the sinu- 
soidal radial-flow path. 

The variation of axial velocities is given in figure 4(c), 
which shows that the axial velocities increase toward the hub 
in all cases, but at different rates. The high value of axial 
velocity at the hub entering the rotor blade allows the use 
of higher turnings at all radii without exceeding the limiting 
value of (Ve,2—Ve,i)/V Ztl or <jC l at the hub. It also helps 


A - 2, approx imafel 


'l %\ 

nS''\ * ' \ for s - r e.) 
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(g) Pressure distributions and pressure rise across rotor. 


to give a more uniform Mach number relative to the rotor 
blade over the blade height. As a result, this type of design 
gives a higher pressure rise and a higher specific mass flow 
than a free-vortex type of design using the same design 
limitations. In order to utilize this advantage fully, the 
variation of axial velocity should be correctly determined. 

The calculation of axial velocity based on simplified i'adial 
equilibrium gives a result close to the zero-aspect-ratio case, 
which is also true in the distribution of other properties in 
this calculation, because in the case of zero aspect ratio, the 
curvature caused by radial motion is negligible and the 
difference in gas properties caused by the radial displacement 
is very small in the nontapered passage. 

The variation of tangential velocities is shown in figure 4(d). 
These velocities in different cases vary in a similar manner 
and the difference of magnitude between them is mainly 
due to the different value of 8 t determined by the different 
values of V Zfi ,JU t in the various cases. 

Figure 4(e) shows the variation of air angles entering the 
rotor and stator blades. The difference between the simplified- 
radial-equilibrium calculation and the case of aspect ratio 
of 2 is significant throughout the whole blade height. The 
simplified-radial-equilibrium calculation gives a value about 
3° lower than the aspect ratio of 2 at the tip of the rotor 
blade and at the hub of the stator blades. 

The variation of Mach number relative to the rotor blades 
is shown in figure 4(f). The simplified-radial-equilibrium 
calculation gives a nearly constant value; whereas the more 
correct calculations show that Mach number actually de- 
creases about 10 percent toward the tip for the case of blade- 
row aspect ratio equal to 2. (This variation, however, is 
only about one-third of that of a similar free-vortex com- 
pressor.) 

The pressure distributions ahead of and behind the rotor 
and the pressure rise across the rotor at different radii are 
shown in figure 4(g). The difference in pressure distribu- 


A - O A ~ 2 A ~ co 



(h) Velocity diagrams at different radii. 


Figure 4. — Concluded. Symmetrical-velocity-diagram and constant-total-enthalpy compressor. 


87359Q — 50 3 
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tions may explain to a certain extent the difference found 
between measurement and the simplified-radial-equilibrium 
calculation. The pressure rise across the rotor is fairly 
uniform in the case of an aspect ratio of 2 and is a desirable 
feature. 

The velocity diagrams at three radii for aspect ratios of 0, 
2, and “ are shown in figure 4(h). If this stage is used as 
the first stage of a compressor, the permissible tip rotor 
speed of this design at standard sea-level conditions is equal 
to 868 and 826 feet per second for A — 0 and A= 2, respec- 
tively. The specific mass flow per unit annulus area, cor- 
rected to standard sea-level conditions, is equal to 41.5 and 
40.0 pounds per square foot per second for ^4=0 and ^4=2, 
respectively. 



GjGt r a ~ Ti 

r t - rn 

(a) Distribution of specific mass flow. (b) Radial displacement across rotor. 



Free-vortex and constant-total- enthalpy compressor. — 

The design constants used are the same as in the previous cal- 
culation. In addition, V& i and V\ 2 are considered to be equal 
to Wej and Ws ,t, respectively, at the mean radius. In this 
type of design, the simplified-radial-equilibrium approxima- 
tion is equivalent to the zero-aspect-ratio case, because, due 
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(d) Variation of air angles. 
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(e) Variation of Mach number relative to rotor blades. 


Figure 5.— Free-vortex and constant-total-enthalpy compressor. 
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to the constant value of rV» in this design, the radial motion 
across the blade does not affect the calculation in the zero- 
aspect-ratio case. That is, the same values of Hi, H 2 , { 1 , f 2 , 
V z 1 , and V z2 occur in both cases and the entire calculation 
is the same. (See also equation (F6).) 

Because the radial motion involved in this type of design 
is mainly due to the compressibility of gas, the difference 
between the zero- and infinite-aspect-ratio cases is not large ; 
hence the calculation for a finite -aspect-ratio case is not 
made. 

The distribution of specific mass flow ahead of and be- 
hind the rotor is presented in figure 5 (a) . Even in the zero- 
aspect-ratio or simplified-radial-equilibrium case with a 



constant axial-velocity distribution, considerable change 
occurs in density, which requires an appreciable amount of 
outward radial motion to obtain the given design conditions 
behind the rotor. Although the amount of this radial motion 
is small (fig. 5(b)), its effect on the variation of gas properties 
is not entirely negligible. Its effect can be seen in the curves 
of figures 5(c) to 5(g), which are somewhat similar to the 
symmetrical-velocity-diagram and constant-total-enthalpy 
design in nature but of smaller magnitudes. 

If this stage is used as the first stage of a compressor, the 
permissible tip rotor speed at standard sea-level conditions is 
equal to 758 feet per second for A= 0. This tip speed is 
about 13 percent lower than that of the corresponding case 
of the previous design. The specific mass flow corrected to 
standard sea-level conditions is equal to 38.6 pounds per 
square foot of annulus area per second for A—0, which is 
7 percent lower than that of the corresponding case of the 
previous design. 

Free-vortex and constant-total-enthalpy turbine. — The 

design constants used in the calculation are : Mi<Cl, U t /a ut — 
0.5, i=0.8, V z ,iJa ut =0A, F 9 , 2 =0, and polytropic 

efficiency at mean radius equal to 0.87. For the simplified- 
radial-equilibrium approximation or zero aspect ratio, 
VJUt is constant and so is V Zt2 /U t , which is found by the 
continuity relation to be equal to 0.877. The same velocity 
at station 2 is used for the infinite aspect ratio, thus making 
the only difference at station 1. The results of the calcula- 
tion are shown in figure 6. 

The distribution of specific mass flow ahead of and behind 
the rotor is shown in figure 6(a). Because of the constant 
axial exit' velocity, the specific mass flow is constant behind 
the rotor. Except for the case of infinite aspect ratio, there 
is an inward radial motion of gas in passing through the rotor 
(fig. 6(b)), the magnitude of which is about two and one-half 
times that in the previous free-vortex compressor (fig. 5(b)). 


A=0 or s.r.e. A-oo 




A t hub 



(g) Velocity diagrams at different radii. 

Figure 5. — Concluded. Free-vortex and constant-total-enthalpy compressor. 



(a) Distribution of specific mass flow. (b) Radial displacement across rotor. 
Figure 6. — Free-vortex and constant-total-enthalpy turbine. 
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(e) Variation of axial velocity ahead of rotor. 



(d) Variation of gas angle entering rotor blades. 



The variation of axial velocity ahead of the rotor is shown 
in figure 6(c). An increasing axial velocity toward the hub 
of about 15 percent would be required for an aspect ratio of 2. 

Figure 6(d) shows the radial variation of gas angles 
entering rotor blades. The difference is only important at 
the hub. In the actual case of an aspect ratio of 2, the 
simplified calculation would give an angle of attack 3° to 4° 
too high at the hub. 

The absolute and relative Mach numbers of gas ahead of 
the rotor are shown in figure 6(e). In the actual case of an 
aspect ratio of 2, the Mach number at the hub is about 
3 percent higher than the simplified calculation. 

Figure 6(f) shows the pressure distribution ahead of the 
rotor. For an aspect ratio of 2, the pressures at the tip and 
at the hub are about 2 percent higher and 3 percent lower 
than the simplified calculation, respectively. 

The velocity diagrams at three radii for the zero and 
infinite aspect ratios are shown in figure 6(g). 



(e) Variation of Mach number ahead of rotor. 


A=0 or s.r.e. A=oo 



At meon rod i us 


*V 2 .K 


(g) At hub 




(f) Pressure distribution. 


(g) Velocity diagrams at different radii. 


Figure 6.— Concluded. Free-vortex and constant-total-enthalpy turbine. 
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SUMMARY OF ANALYSIS AND CALCULATIONS 

In axial-flow turbomachines, radial motion of gas occurs 
because of the tapering of the passage walls and the variation 
of gas conditions across blade rows specified in the design. 
The direction and the magnitude of this radial flow depend 
on the type of design, the tapering of the passage walls, the 
blade-row aspect ratio, and the Mach number of gas flow. 
Even in the free-vortex type of design employing nontapered 
passage walls and requiring no change in velocity distribu- 
tions from stage to stage, an appreciable amount of oscillatory 
radial motion occurs within the stage. 

This radial motion gives an additional term to the ordinary 
radial-equilibrium equation. In the free space between blade 
rows, this additional term is approximately equal to the 
product of the square of axial velocity and the curvature 
caused by the radial flow. Depending on whether the 
curvature is positive or negative, the radial pressure gradient 
caused by the whirling motion of gas is decreased or increased, 
respectively, by this additional term. 

The determination of this radial-flow path requires a long 
process of step-by-step calculation. It is found, however, 
that a sinusoidal radial-flow path gives an effect on the 
radial variation of gas condition between blade rows as 
small as possible without discontinuity in the curvature of 
the streamline. It may therefore be desirable to prescribe 
this simple radial-flow path in the design. Also, inasmuch 
as it represents the major harmonic of the radial-flow path 
that may exist in any design in which the blade loading is 
relatively uniform in the axial direction, the calculation 
based on this simple radial-flow path gives good, approximate 
results. 

Methods of solution for the limiting cases of zero and 
infinite blade aspect ratios and a simple approximate solution 
of the radial displacement across the blade row having a 
finite aspect ratio are also obtained. 

The analysis made of the maximum compatible number 
of the degrees of freedom in specifying the radial variations 
of gas properties in stations between successive blade rows 
of a turbomachine shows that the designer is free to specify 


one such variation at each of the stations. The various 
ways to use up these degrees of freedom and the resultant 
types of design obtained are discussed. 

The usual method of calculation, which neglects the radial 
motion, gives results close only to the case in which the axial 
length of the blade row is much larger than its radial length, 
and is not good for the case of a finite blade-row aspect ratio. 
The difference between the results obtained by the usual 
method and the method suggested herein is found to be 
quite large in a design employing constant total enthalpy 
and a symmetrical velocity diagram along the radii. Calcula- 
tion made for this type of compressor, using the same limit- 
ing Mach number, same limiting turning, same axial velocity 
at the mean radius, and for a blade-row aspect ratio of 2, 
gives the following differences between the usual and the 
suggested method: 

1. The radial variation of axial velocity ahead of the rotor 
is 13 percent for the usual method and 28 percent for the 
suggested method, and the radial variation of axial velocity 
behind the rotor is 53 percent for the usual method and 
40 percent for the suggested method (all expressed in terms 
of their values at the mean radius). 

2. The air angles differ from 1° to 3° at the hub and at 
the tip. 

3. The radial variation of Mach number relative to the 
rotor blade in the usual method is 9 percent lower than 
that in the suggested method. 

4. The radial variation in static-pressure rise across the 
rotor is 13 percent for the usual method and only 2 percent 
for the suggested method. 

5. The mass flow in the usual method is 4 percent higher 
than that in the suggested method. 

6. The allowable rotor speed in the usual method is 
5 percent higher than that in the suggested method. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, January 1, 1949. 





APPENDIX A 

SYMBOLS 


The following symbols are used in this report: 

A aspect ratio of blade row, r * j / h 

velocity of sound 
lift coefficient 

specific heat of gas at constant pressure 
specific heat of gas at constant volume 

differentiation with respect to time following motion 
of gas 

external force per unit mass of gas (general case); 
circumferentially uniformly distributed blade force 
per unit mass of gas_ (axially symmetric case) 
radial component of F 
axial component of F __ 
tangential component of F 
form of radial-flow path 

mass flow per unit flow area perpendicular to axis 
of turbomachine 

form of radial-displacement distribution 

V 2 

total enthalpy per unit of gas, h -\ — 

7 ) 

enthalpy per unit mass of gas, w-f— 
constant 

thermal conductivity of gas 
axial length of blade row (fig. 3) 

Mttch number of gas 
m mass of gas 

n polytropic exponent of actual expansion or compres- 
sion process of gas 
p static pressure 

Q heat input to unit mass of gas along its path of motion 
per unit time from neighboring gas particles in 
general case or from blade and passage wall in 
axially symmetric case 
R gas constant 

r radial distance measured from axis of turbomachine 

r 0 mean radius of radial-flow path (fig. 3(a)) 


s entropy per unit mass of gas 
s. r. e. simplified-radial-equilibrium approximation 
T absolute stream temperature of gas 

t time 

U magnitude of U 

U vector velocity of blade at radius r 
u internal energy per unit mass of gas with 0° absolute as 
base temperature 
V magnitude of "V 


V 


V r 

v, 

V e 

W 

W 

W e 

y 

z 

a 


y 

A 

5 


absolute vector velocity of gas 
radial component of 'V 
axial component of V 
tangential component of V 
magnitude of W 

vector velocity of gas relative to rotor blade 
tangential component of W 

maximum radial displacement over blade height 
distance along axis of turbomachine 
angle between absolute velocity of gas and axis of 
turbomachine 

angle between relative velocity of gas and axis of 
turbomachine 
ratio of specific heats, c P /c v 
radial displacement across rotor, r 2 —r 1 

y y 

dimensionless turning , i2 r / 1 

r \Ul ,t 


f angular momentum about 2 -axis per unit mass of 
gas, rV e 

ij small-stage or polytropic efficiency 
8 angular coordinate measured from some fixed radial 
line 

/z absolute viscosity of gas 

p mass density of gas 

<r blade solidity 

$ dissipation of energy due to viscosity per unit volume 
of gas per unit time 
c p function 

X angular coordinate measured relative to rotor 
w angular velocity of blade 

Subscripts : 

1 ahead of rotor 

2 behind rotor and ahead of stator 

3 behind stator and ahead of next rotor 

c satisfying continuity equation 

e satisfying radial-equilibrium and total-enthalpy equa- 

tions 

h hub 

i any station between two blade rows 

j any station 

k station short distance downstream of station j 
l limiting value 

m at mean radius 


n used with r to indicate radius where maximum radial 
displacement occurs 


s simplified-radial-equilibrium approximation 

t tip 



APPENDIX B 

DERIVATION OF EQUATIONS 


From equations (4), (7), and (9), and the relations 


R—Cp—c t 

and 

(Bl) 

Cj> 

7=— 

Cr 

(B2) 

there is obtained 


dh= d( 2 A 

7—1 \p/ 

(B3) 

From equations (8) and (B3), 


ffi-Fi'* (?)+■* G) 

(B4) 

From equations (9) and (10), 


Tds=dh~— 

p 

(B5) 

By use of equation (B3), 


T< h=-^rd(e)-& 

7—1 \p/ p 

(B6) 

By use of equation (7), 


d (s\ 1 ip y dp 

\RJ 7—1 P 7—1 P 

(B7) 

=^r x d(loz,f) 

(B8) 

d (log, (log, p) 

(B9) 

~ y —l ^( lo ^ T) d (log, p) 

(BIO) 


Equation (1) may be written as 

V-F+j \ log ‘ p=z ° 

Combining with equation (BlO) yields the following form of 
the continuity relation: 


From equation (2) and the following relations 


(la) 


and 


D ^=^+(V.V)V 


(V-V)T=5 VF‘- FX{VX V) 


there is obtained 


dF 

dt 


™ 1 


-FX (VXV)+| VF 2 =F~ [v 2 F+§ V(V-F)]+ 
\ {2[(V m ).V]T+(V m )X(VXF)-| (V-FHVp) j (Bll). 


1 „ T* 


From equations (B4) and (B6), 

1 VF 2 +- Vp=VH 

2 p 7—1 


+ivj. 


=VH- TVs 

Combining with equation (Bll) yields 

vH—F+ tvs+vx (vxf ) -^+- p [v 2 f+| V(V.?)]+ 
~ p | 2 [(Vp)-V] F+(Vp)X (VX V)-\ (V-F) (Vp) | 
From equations (8) and (B5) , 




3 V(V-V)J 


+ 


PJP rp IDs , 1^ Dp ,'Tt DV 

Dt ~~ 1 Dr P Dt^ v ‘ Dt 

= T m+\t+ T ^(h v t +I §) 

By combining with equation (2), 

w= r §H! +T? -( F +?[ v: 

\ { 2 [(Vp) -V] F+ (Vp) X (VXD) -| (V-F) (Vp) j) 
When combined with equations (3) and (10), 

=«+*■ ^+ V - (*+? [v^+| V(V.F)]+ 

~ | 2 [(Vp)-V] F+ (Vp) X (VX F) — | (V*F) (Vp)j)(3a) 


DU 

Dt 


From equations (3) and (10), 

Ds 
Dt 


Zm=<2+f 


Dt~T' u p 


(3b) 


For steady axially symmetric flow of nonviscous fluid, 
equation (3 a) reduces to 


From equation (11), 


^Q+FV 


F-V—F-U 
= FeU 
= rFou 


(B12) 


(B13) 
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From equation (15), 


rF.= V,*^ + V, 

D(rV„ 


Dt 


(B14) 


Combining equations (B 12 ), (B13), and (B14) yields 


DH_ n , D(rV e ) 

Dt~ Q+a ~DT 


(18) 


For steady axially symmetric flow of a viscous fluid, 
equation (18) is obtained by applying motion and energy 
equations to a mass system with a fixed control surface, 
as shown by the solid lines in figure 2 . Under steady 
axially symmetric flow, the mass inflow drrij in time dt is 
equal to the mass outflow dm k in time dt, the state of gas 
within the control surface is unchanged, and the state of 
gas at stations j and k is constant with respect to 9. By 
equation (Bll), the sum of the tangential blade force and 
the tangential viscous gas forces exerted by the surrounding 
gas particles on the system is equal to 


1 D(rVe) 
r Dt 


dm i 


The torque about the 2 -axis exerted by these tangential 

forces on mass dmj is therefore simply dm, and the 

work input to mass dmj by these tangential forces in time 
dt is equal to 

^ dmj u dt = u[ ( rV a ) * — (rV g ) } \dm 3 

In passing from station j to station k , in addition to receiving 
this work input, the gas particles are doing work against the 
axial and radial viscous forces exerted by the surrounding 
gas particles. This negative work is usually small, however; 
if it is assumed that the heat generated by the frictional 
work is added back to the same gas stream, the heat addition 
cancels the negative work, and the energy equation for 
steady flow gives 

(H k — Hi) dm, — k Q dm } dt4-u[(rV e ) k — (rV 9 ) } ]dm } 


DH n , D(rVg) 

-Di=i + “—Dr 


(18) 


where Q denotes the rate per unit mass at which the gas 
stream sheet is receiving heat from external source through 
blades or other passage walls. 

When equation (19) is given, the entropy change can be 
obtained in the following manner: From equation (BIO), 


D 

Dt 


(r)~- 


1 Dt 


D T D . 
Og* T— w, lo£ 


Dt 


(B15) 


But by equations (19) and (7), 


D . 1 D , T 


Substituting into equation (B15) gives 

DS r> ( 1 1 \ D y rp 

Dt~ R \y - 1 n-l)Dt logeT 


=R 


71 y -j rji 

(71-1) (7-1) Dt 10ge J 


D 


(B16) 


For steady axially symmetric flow, equation (B16) reduces 
to 


Ds „ n — 7 
Dt~ U (71-1) ( 7 


^ I y(r,| io g.r + F,|iog < r) m 


For successive axial stations j and k a short distance apart, 
equation (Bl 6 ) gives 

s k (r k )—s j (r j )=R ~i j t loge T k (r k )—\og e Tj(r } )] 


=R 


71 — 7 


(71— 1) ( 7 —I) 


log, 


T k (r k ) 

Tj(rj) 


(B17) 


Inasmuch as the temperature change between the two 
successive stations is small, if the enthalpy is measured with 
0 ° absolute as the base temperature, the temperature ratio 
can be considered equal to the enthalpy ratio: 


_ 

Tj hj, 



(B18) 


Substituting equation (B18) into equation (B17) gives 

7 2 

TT ’ k 
til- 


Sk(r k ) -s } (r } ) =R 


71—7 


(71—1) (7—1) 


log 


_Z* 2 

K cy 

TT v J 
2 


( 20 a) 


The density ratio between the two stations is obtained 
from equation (BIO) 


Pk_ 

Pi 


Combining with equation (B18) yields 


(B19) 


Pi 


H k - 


V k 2 


H- 


V/ 


' R 


(B20) 


Substituting into equation (22) gives 
1 


V 2 , k (n k - V -^y 1 e R r k dr k ~V z , } ' e * r } d rj 

( 22 a) 



APPENDIX C 

DETERMINATION OF RADIAL DISPLACEMENT BY USE OF CONTINUITY EQUATION 


Equation (22) may be written as a linear differential equa- 
tion for r 2 as a function of r x , provided G 2 (t\) is known: 


dir, 2 ) 




r 2 


dr x u Giin) 
* G x jr x ) 


J r i G 2 (r x ) 

when divided by r 2 , t 2 , 


2 r x dr x +r 2 . t 2 


(Cl) 


r~Y=i— Y f r 

V2,l/ V2,(/ J r 1 


r *.‘ G x (r x ) 2 t\ dri 
G, (r x ) r XJ r ut 


(C2) 


If G x and G, are unknown and only G x jG lit and G 2 /G 2 , t are 
known, a modification is necessary : 


/ >A 2 — 1 _/ rjSS Gut C T ut 

V2 ,t) V2,(/ G 2t t Jr l 


G x , , 

G 2 , , 

7T W 

1*2,1 


r i dr x 

r7,t r\i 


(C3) 


The value of G\ t tlG 2tt is found by the condition that total 
mass flow at stations 1 and 2 is the same: 


Hence 


fe) 


2 Gu,. 


/rjA’ 

Va.i/ 


g 2 


r 

J r u> 


G± , , 

.t Gut 2 A 

IrA 


r 


G x 


(r.) 




_ _ A dr x 
G 2 {r yutrx,t 




,1 Out r ± dr, 

1*2 A 


(C4) 


(Co) 


APPENDIX D 

EQUATIONS FOR SIMPLIFIED-RADIAL-EQUILIBRIUM CALCULATIONS 

This equation holds for all stations, provided the appropri- 
ate values of (Ve.t/dt) are used. It follows from equa- 
tion (D2) that at each station 


Equations to calculate distributions of gas properties at 
three stations of a typical stage under the simplified-radial- 
equilibrium approximation for a few types of design are 
given. 


GROUP I 

Free vortex. — For this design, 

dU 


dr i 


=0 


(Dl) 


When the inlet total enthalpy is constant with respect to r 
and the radial variation of entropy is negligible, equation (14g) 
reduces to 

dV Zii 


dr t 


'-0 


(D2) 


The variation in tangential velocity is, by equation (Dl), 


Ve.t — V 9, i.t 


r { 


(D3) 


At each station, by using equations (30) and (D3), 


7 |/? s 

r pt y dr 


OAaL 2 

9,3 


When the preceding equation is integrated from r to r„ an. 
the relation 

frp 

p 

is used, there is obtained 


°=V- 

ri-^eyRw-r 


(D4) 


G_p 
G t pi 


(D5) 


The radial position of gas at station 2 or 3 can be obtained 
by numerically integrating equation (C5) using distributions 
of specific mass flow given by equation (D5). An alternate 
method is to expand the right side of equation (D4) into a 

binominal series. Because 7 9 1 — * J * s usua ^Y 

less than 0.15, three terms will be sufficient. Let p r repre- 
sent the average density in the annulus between r and r„ 
then 


- P 

J r 


2 irpr dr 


Pt x (r, 2 — r 2 )p ( 


2-7 

4 


/F fl ,An lQge O | 2 7 
UJJ i-0 2 


(«)’ 


(D6) 


and inasmuch as 

Pr, 1 V,.i (r ut 2 -r x 2 ) =-p rA V z . 2 (r 2 .: 2 -r 2 2 ) 
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Pr,l 1 . ? — ^l 2 ) V 2,2 Pft.l — 7*1, ft 2 ) 

Pr .2 (r 2 , t 2 —r 2 2 ) ~V z ,i~~p hA {r 2 J—r 2 , ft 2 ) 



The change of total enthalpy across the rotor is 
H 2 Hi=w(^2 — ft) 

Then 

. , w(f 2 — fl) 

Hr 1+ h x 

For compressor, 

TT 1 7 

H 2 , , r lth r,.» Ve ^ U ht 2 

Hi - 1 + r ut U UI V z ,i,h Hi 


where 


^ Ve^-Ve.un 


the quantity — — wy~~ is t° be chosen by the designer. 

V z.i.ft 

Inasmuch as 

Hi=h ut +^f 


(D7) 


and 


c h 


=^Ti+\ (V zX f+V e ,i,i 2 ) 

2 _/fu V 

\M l ) 


=^2 W*f+(U l . t *-v, A .d*] 


where M, is the limiting Mach number to be chosen by the 
designer, 


H* 

Hi 


= 1 + 


- 2 - fl t 7 V 

r V r y 9 , 2 .h—V 6 ,l,h 
7 l,ft v z,l,ft '1 .ft 

r l,t Hut Fg.l.ft 


i r/ ^.i.i V , (, , i VfVz.iA^/Ve.iAn 

( T - 1)M, 2 LV U ut ) + V U ut ) J+2 |A U i.t ) + VIW J 


(D8) 


The pressure distribution at each station is obtained by 
raising its density distribution (equation (D4)) to the power 7 . 
The pressure changes between the stations at different 
radii are obtained by combining these pressure distributions 
with the pressure change across the rotor at the radius where 
the value of the polytropic exponent is known or assumed. 
The angle that the gas velocity makes with the axis of the 
machine at any radius is obtained from the known tangential 
and axial velocities. 

Symmetrical velocity diagram. — For the nontapered pas- 
sage, r x ~r 2 =r and from equations (39) and (42a) 


dU 

dr 


d^2 

dr 


j. wr 2 r ht 8,U ut 
2 2 




(D9) 


where the minus sign is used for station 1 and the plus sign 
for station 2. Inetgration from r h to r yields 


={ 


1 + 


7-1 U{ 

8 a h 2 






(Dll) 


where the minus sign is used for station 1 and the plus sign 
for station 2. 

For the case where the inlet total enthalpy is constant with 
respect to radius and the radial variation of' entropy is 
negligible, the variation of axial velocity is obtained from 
equation (14g): 


J 

Ve,i _ 1 _r tnrij 

U ut 2 r us 2 r 

Zo = l r_ 

U ut 2 r lit ' 2 r 


(DIO) 


When equation (DIO) is substituted into equation (30), 


7 



dp__ai 2 r _ 
c>r 4 


Sj ^-+\ MU,.,)’ 


1 

y»3 


v t- 

2 dr r 2 ^ dr 


Integration from r h to r gives 


(D12) 


(D13) 


where the plus sign in the last term is used for station 1 and 
the minus sign for station 2. 



APPLICATION OF RADIAL-EQUILIBRIUM CONDITION 

When both Ve and V z are known, the radial variation of 
density can also be obtained by applying equation (B20) at 
the station: 

i 


p__ 

Ph 


H (Vg 2 +V z 2 ) 
H-\ (V^+V;,* 2 ) 


7-1 


(D14) 


In the tapered passage, the gas is assumed to flow in 
conical surfaces, which gives the value of r 2 as a function of 
r x . Equations (39) and (43) give the distributions of and 
Ve as shown by equations (E5) and (E6), respectively, given 
in appendix E. The distribution of axial velocity at station 2 
is the same as that given by equation (E8). The density 
distributions can be obtained from equation (D14). After 
these distributions are known, the distribution of specific 
mass flow G is known and the radial displacement is found 
by using equation (C5). 

In compressors of this design, the maximum value of 

— j s usually at the hub. Its value there is to be 

set by the designer. Then 
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Integrating from hub to radius r yields 

Vz. 2 =V z ,i,h 2 —2 K 2 {r l 2 —r Xth 2 ) 


or 

With radially constant work input to the rotor, and a non- 
tapered passage, 


5*= 


s- s- Ve, 2 ,h— Ve, l.ft T r 

' h ~r uh U ut - Vz.i.n u ut 


and 


r ut 


?*i ,* V z , i,ft r lt ft 


r ^Ve.uH~Ve,un 


n, t Ui.t 




(1)15) 


In this type of design, the limiting Mach number is usually 
at the hub. Hence the denominator of the last term of 
equation (D8) should be replaced by 

WS-WMSWSI 


(y- 

The rest of the calculation is the same as in the previous 
design. 

Wheel-type tangential velocity in front of rotor. — When 
the case of constant total enthalpy is again considered at 
the inlet and the radial variation in entropy is neglected, 
with 

V,, 1 =K l r l =V a .utP- (D16) 


r U i 


equation (14g) gives 


dVzl 

dr x 


= — 4 K x r x 


and 


r 2 F 9 ,2=^Lir 1 2 4-5 { r lit C7 1 ,j (D18) 

V,.,., 


4^J 1 log.i 

u x , t r l,h 


(D19) 


Equations (D17) and (Dl8) show that the axial velocity 
rapidly decreases with radius at stations 1 and 2. If K x in 
equation (D16) is chosen to be to/2, the difference between 
this type of design and the previous one is very small. 


GROUP II 


Untwisted rotor blade. — Equation (52) gives 
Vg i — C or i 




s =tan &i=Ki 


(D20) 


where i=l, 2. When equations (14g) and (D20) are used, 
the following relation is obtained: 

K ,, orndVe.t , Vg,i 2 <aV 9 ,i dHt „ dSj cSri 

dr t + r, K 2 ~ dr i li dr i K? 

(D21) 


or 


[(14 -Ki 2 )Vz,i+KiOiri] VzJ+ZKtuVz. 


, i -T i -p~2u 2 r t 
dv i dv i 


(D22) 


Either equation (D21) or (D22) can be solved by a standard 
method of numerical integration. Equation (D20) is then 
used to find the remaining velocity component. In equations 
(D21) and (D22), dH>Jdri in later stations, except at the 
station ahead of first rotor, is, in general, not equal to zero 
even if it is equal to zero at the inlet. These ratios are to 
be determined by using equation (18b). The term contain- 
ing entropy in the equations may be significant in the case of 
cooled turbine. 



APPENDIX E 

EQUATIONS FOR ZERO- AND INFINITE-ASPECT-RATIO CALCULATIONS 


The method of calculations is given for the two types of 
design used in the numerical examples. 

FREE VORTEX 

Zero aspect ratio. — In this type of design, the zero-aspect- 
ratio case is the same as that of the simplified-radial- 
equilibrium approximation. 

Infinite aspect ratio. — By equations (31) and (B20), 


V:,2 

Vz,i 



2H l -(Ve, 1 2 +V lA 2 ) 
2H 2 — (T 7 0,2 2 + Vz'i) 


(El) 


An additional relation between V z ,\ and V z , 2 is necessary 
in order to solve the equation. In the section Limiting case 
of infinite aspect ratio, two equations are suggested. For this 
design, equation (33) gives 


or 


dV zA | dV z . 2 n 
dr ^ dr -u 

Vz. i + Yu = constant 


and equation (34) gives 


or 


Vz. 


dV z . 

dr 


l + V 2 , : 


dV z 


dr 


V z .i 2J r Vz.i~ constant 


(E2) 


(E3) 


Also from equations (F4) and (F5), when the square term 
in A„ is neglected 

F** 2 t i F*" 2,2 constant (E4) 


(4) Plot V 2 ,j and V z , 2 against rjr t , and obtain V Zjl and V Zt2 
at the values of r/r, desired. 

When the distribution of axial velocity is known, the den- 
sity variation at any station is obtained from equation (D14). 
The pressure variation at each station is obtained by raising 
the density ratio to the y power. The pressure changes 
across the stage at different radii and the air angles are ob- 
tained in the same manner as in the simplified-radial- 
equilibrium calculation. 


SYMMETRICAL VELOCITY DIAGRAM 


Zero aspect ratio. — With radial displacement not equal to 
zero, the equations for tangential velocities are different 
from the expressions of equation (DIO) . From equations (39) 
and (43) 


dr x dri 03 1 


and 


fi 


£2 


2 2 ~ 

wr t 2 r ut S t U ut 

2 ' 2 


V$, 1 Ti r i 

U ut r\ Ui.t 2ti,( 2rj 

V 8 , 2 _ £2 __/ r i 1 r ut S, \ r 1 

U\,t r 2 Ui.t \2r 1>{ 2 iq / r 2j 


(E5) 


(E6) 


From equation (14g), neglecting the radial variation in 
entropy, 


t / dV z .j 1_ .. d£j 

2,1 dr t ~ r t 2 dr { 


For station 1, from equation (E5) 


The three preceding equations give practically the same 
results. 

A convenient procedure of calculation is as follows: 

(1) In order to compare the result with other oases, the 
same value of V zA . m may be used. From equation (El), 
Vz t2 ,m is determined. 

(2) Insert these values in equation (E2), (E3), or (E4) to 
obtain the constant in the equation. 

(3) Assume a number of values of F 2>1 ; obtain V z , 2 by the 
same equation. Then use the following equation, which is 
obtained from equations (El) and (D3), to solve for rjr t : 


(Ve.2 ,,) 2 ~ 

/ *2 *1 
1 e R 

v z ~) 

,T-I 

I ( V e , M ) 2 

2 H, 

1 

1 1 


’■‘>1 

f *2“' 

r 

~V zA } 

VzJ 

€ 

J 

7 


xr dV z . 1 1 / wrr ri TT 2 /l ^ 1 \ 

~di\ = ~r? VT 2 ) “ r ' = - U '-‘ \2 rj 

Integrating from r h to r l yields 

fVz A \ 2 _/Vz A . h \ 2 fi 1 (r, 2 r,,A 

KwJ -\TT 7 ,) +5,Iog ' (E7) 

which is the same as equation (D13). For station 2, 


dVz 

dr 


,2 1 / cor i 2 | r U1 5,U 1 

~ r 2 2 \ 2 + 2 ) 


wri 


dr 1 
dr 2 


Vz,2dVz,2=~~2 (c W + SM.fldn 


2 r 2 


(B 
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Integrating from r ltfl to r x yields 


(t + x) d (B 

n.t \ r x ,J 


(E8) 


which differs from equation (D13). 

The density distribution at station 1 is the same as the 
simplified radial approximation, whereas that at station 2 is 
obtained by using equation (Dl4). The solution of this 
case is a process of successive approximations. Values of 
r 2 (ci) obtained in the simplified radial approximation can be 
used here as the starting values. Then the distributions of 
Ffl.o, F 2i2 , P 2 , and f? 2 are calculated from the preceding equa- 
tions, and new values of r 2 (/*i) are computed from equation(Co). 
Usually, only two or three cycles are necessary to obtain 
the correct value, because the difference between this case 
and the simplified-radial-equilibrium approximation of this 
type of design is small. 

Infinite aspect ratio. — The first equation for the condition 
Gi—G 2 is the same as equation (El). The second equation 
necessary to solve this case is a little more complicated than 

that in the previous type of design because ^-5^0- If equa- 


tion (34a) is used, 


T T dV zA | T7 - d\ Zl 2 
Vz - l ~dT +Vz ' 2 


dr 


__I ( y- 

r 2 \ dr ' dr J 


Then 

or 


— ~~ (U + fo) — — urr 
V z , i 2 + F 2 ,2 2 = — w 2 r 2 + constant 

Vi, l 2 , V 1,2 

U t 2 + uf 


(E10) 


In order to compare the result of this case with other cases, the 
same value of V ZAm may be used. Then from equation (El), 
V z , 2 , m is found, and the constant in equation (E10) is evalu- 
ated by using this set of V Z Lm and V z2 ,m- A few values 
of V ZA are assumed at any other given radius, with corre- 
sponding values of V z 2 obtained from equation (ElO). The 
correct values of V ZA and V Zt2 that will satisfy equation (El) 
are obtained by interpolation. 

After the distribution of axial velocity is known, the den- 
sity distributions are obtained from equations (Dll) and 
(D14), and pressure distributions, total enthalpy change, 
and air angles are obtained in the same manner as before. 


APPENDIX F 

APPROXIMATE VALUE OF RADIAL DISPLACEMENT ACROSS BLADE ROW HAVING FINITE ASPECT RATIO FOR GENERAL CASE 
IN WHICH Ti AND f 2 ARE PRESCRIBED IN DESIGN AS FUNCTIONS OF n 


In the latter part of appendix E, distribution of axial veloc- 
ity is expressed in terms of known H, r h and r 2 {r x ). Alter- 
natively, this distribution can be expressed in terms of radial 
displacement and its value determined by the simplified- 
radial-equilibrium calculation, for which A e =0. For a 
nontapered passage, it is seen from equation (14c) that 


I/- dV z , i Tr dV z ,i, s 

v 2 ,1 v 2 ,l, s 


+ 


i ay 


V 2 ,1 2 


(FI) 


and 

T . dV 

I 2.2 


dr 


dr i z,i ' s dr y 

z , 2 T7 dV 2 , 2 ,s A c / 7r\ . g dr 2 , . d ^2 ( 1 1 \ 

t * ,2,s dry 2 XL/ ^ z ' 2 ' dr i ^ 2 dry y"i 2 r i) 


(F2) 


By substituting (r 1 + A e ) for r 2 , expanding ^1 -j-^A in a bi- 


nomial series, and neglecting terms of greater order than 
, equation (F2) becomes 


V z 


dV z , 

dry 


'=F, 




+ 


If A e (?‘i) is known, equations (Fl) and (F3) may be solved 
as linear first-order differential equations in V 2 ,i 2 and V z /, 
respectively, giving (omitting the subscript 1 on r) 


V z d=e 




dr + V z ,y,„ 
(F4) 


and 


where 


| £ [£ V z . 2 , 2 +^ £ f 2 2 -^ £ dr+V z , 2 ,J | 


(F5) 


*>(U=£ A e dr+\ {(A e ) 2 -[A e (r ra )] 2 } 

and subscript m may here refer to any radius between hub 
and tip. 


For the limiting case of zero aspect ratio, the last term in 
equation (Fl) approaches zero so that V zx =V zls , whereas 
V, 2 is obtained by integrating equation (F3) with the third 
term neglected. 
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J,' £ (J tf) dr- 



(F6) 


where VJJj is a mean value of V z , ] . s 2 between r and r m , the 
mean depending on the choice of the function g(r) . If the ap- 
proximation is made in letting V Zi2 , m = V Zl2 ,m,s, equation (F5) 
may be written as 


When equation (F4) is integrated by parts and A e (r t ) is 
replaced by y e g(r) as in equation (35), there is obtained 




Si’- 


<pi « 


OW-V., 


si 


Ve<p |W e* r 


in which 


Vi (r)=Vj(r)dr 


-Si 


y e <Pi (r) 


dr 


If it is desired to compare the general case with other cases on 
the basis of the same V 2% hm , then V 2iUm = V ZtUmiS . By the use 
of the mean-value theorem of integral calculus, the preceding 
equation can be written as 



V z ,i-V z . 2 J=V z , 2 J [e (I) 2 ^ (r > -i]+y e<P3 (r)-yM(r) 


(F8) 


where 


Mr)=^p-=Mr)+k y e {[0(F)] 2 — [g(r m )] 2 } 

Ua & 


Mf=e 


and 

<Pi(r)=e 



The change in the distributions of V z ,i and V z%2 with the 
maximum displacement y e for a given g(r) is now determined 
by differentiating equations (F7) and (F8), assuming that 
r m , V z .i, s 2 , and V z . 2 , s 2 are independent of y e 


*w los ‘ v °-'- 


TW/ttV 

: \l) 


(r)e 


(- 1 )'- 


•^i( r ) 


(F9) 


2 4; 

By subtracting equation (FlO) from (F9) and neglecting three small terms containing y e /V z 2 


<P* (r) y 2 d<p 4 (r) 


V, 2 2 V z . 


dy 6 


d . V z ,: 1/ttV ,Jv z , Us 2 G)^ l(r) 
W ] ° St V^ = 2 W ^ l(r) \VM + 


V: 

v z 


z?-®' 


VeVlV) 


1+2/. 


[g (r)} 2 — [g(r m )f ) I 1 <P3 (r) 


Mr) 


2 V z 


(FlO) 


(Fll) 


The equation of continuity, equation (22), may be written 
as 

Sfe-[ 1 +*^ < F12 > 

by replacing r 2 with ri-\-y c g(r). It is here assumed that the 
displacement A c for the continuity equation has the same 
form as, but may differ in magnitude from, A e . If the 
variation of pi/p 2 (r) with y c is neglected, that is, the density 
distribution is assumed to be determined primarily by the 


tangential velocity distribution, differentiation with respect 
to y c for a given g(r) gives 


A 

dy c 


i0g, tT~ = 
v z ,2 


gif) 
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-+- 


Tr a{r) 


1 -fVc 9 -— 1 +ycfc9(r) 


(F13) 


If the same distribution of V zA and V z , 2 satisfies both the 
continuity and equilibrium equations, y c is a function of y e 
determined by the differential equation, which is obtained 
by dividing equation (Fll) by equation (Fl3): 


W = \ 

dy e 2 g(r) 
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(0 


* l(r) [vZ7 W /~G) ^ (r) j 1+ye MfAzMAl 
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(F14) 


i+2/c 9 -^r i+y c ^[g(r)] 
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In order to evaluate equation (Fl4), the form of g(r), 
which is implicit in the equations used, must be found. An 
order-of-magnitude result may be obtained, however, by 
equating the right side of equation (Fl4) to a constant — K 2 
and determining the value of the constant from the boundary 
conditions on g{r). Because this assumption involves set- 
ting dy c /dy e equal to a constant, it is equivalent to the 
assumption already stated that for the selected A e =y e g(r), 
the corresponding A c differs only in amplitude. In order to 
obtain the order-of-magnitude result, the right side of equa- 
tion (F14) is simplified by 

(a) Setting the first two terms in the bracket equal to 2 
, (b) Considering the terms involving y c and y e negligible 
when they are compared with unity 
(c) Ignoring the last term in the bracket because it 


contains V Zt2 2 \j^j in the denominator (If equa- 
tion (F14) is written in terms of rjr t instead of r, 
the term becomes > which is about 


250 for A— 2 and —= 0 . 6 .) 

r t ' 

As a result of this simplification, equation (F14) becomes 


iy< (l) »' (r) 

dv ‘ £,w 


-IC- 


(Flo) 


Rewriting equation (F15) gives 


<0 


K 2 




(F16) 


When equation (F16) is differentiated with respect to r and 
the relation ^ <piir)=g{r) is used, 


first eigenvalue of this boundary-value problem must be 
taken. A satisfactory approximation to this solution may 
be, obtained without involving Bessel functions by replacing 
g( r )/r in equation (F16) by g{r)jr m , differentiating, and 
solving: 

r 

g(r) —e 2tm (Ki cos lr-\-K 2 sin Ir) 

where 

' = V(z) 

The boundary conditions are determined by using the first 
eigenvalue for K, 


g(T)=K$e 2r " sill -7 — ~n (FI 7) 

r t o» 

and therefore, 

/7rV 1 _ T 2 1 7T 2 

\LJ K 2 (r,— r fl ) 2 ' 4r m 2 (r t — r h ) 2 


(This approximate equality is correct within 1 percent for 
(Th/r t )^> 0.5.) Substituting this result in -equation (Fl 5 ) 
gives 


dye Tv-2 

dy e ~ 



(F 18 ) 


In this very rough approximation, dyjdy e is therefore equal 
to minus the square of the aspect ratio. By integrating 
equation (Fl 8 ) and by letting y CiS equal the value of y c 
corresponding to y e =0 (simplified-radial-equilibrium approx- 
imation), ; 

yc=y c . s ~ A 2 y e (Fid) 


dr 2 g(r) + r dr g ^ + \(lk} r 2 J 9 ^ ~~ 0 


A solution corresponds to y c —y e =y y which when substituted 
into equation (F18) gives ’ 


This equation gives g(r) as a Bessel function of the first 
order and argument (7 r/LK). The value of (tt/LK), and 
thus of K, is determined by the boundary conditions g{r h ) = 
g(r t )—0. In order that g(f) have a single maximum, the 


Ve 


A— 


l-M 2 

A r , 


l+A 2 


(37) 



APPENDIX G 


EQUATIONS FOR FINITE-ASPECT-RATIO CALCULATION 


In the numerical example of the symmetrical- velocity- 
diagram and constant-total-enthalpy compressor, computa- 
tion is made for a blade-row aspect ratio of 2, with a pre- 
scribed simple sinusodial radial-flow path. Inasmuch as 
the term containing radial variation of entropy is not con- 
sidered, substituting equation (E5) into (14c) gives 


and 


<e ) 2 



(Gl) 


dV z , 2 2 
dr^ 


(£)V,- 



From the relation 


r t — r 2 _ 1 / A XV- 1 — ^ 
L 2 ^LrilVi r ‘ / 


and integrating equations (Gl) and (G2) from r n to r, there 
results 


and 



(G3) 


(G4) 


With this set of equations replacing (E7) and (E8), the 
rest of the calculation is the same as in the zero-aspect- 
ratio case. Equation (37) is used as starting value and by 
appropriate interpolation after two cycles of calculation the 
third or fourth cycle usually gives sufficient accuracy. 
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